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The presence and evolution of apparent horizons in a two-parameter family of spherically symmet- 
ric, time- dependent solutions of Brans-Dicke gravity are analyzed. These solutions were introduced 
to model space- and time-varying gravitational couplings and are supposed to represent central 
objects embedded in a spatially flat universe. We find that the solutions possess multiple evolving 
apparent horizons, both black hole horizons covering a central singularity and cosmological ones. 
It is not uncommon for two of these horizons to merge, leaving behind a naked singularity covered 
only by a cosmological horizon. Two characteristic limits are also explicitly worked out: the limit 
where the theory reduces to general relativity and the limit where the solutions become static. The 
physical relevance of this family of solutions is discussed. 
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I. INTRODUCTION 

Varying "constants" of nature, first hypothesized by 
Dirac [1], can be implemented naturally in the context 
of scalar-tensor gravity, in which the gravitational cou- 
pling becomes a function of the spacetime point [2, 3]. 
String theories [4] contain a dilaton field coupling non- 
minimally to gravity which mimics a Brans-Dickc-likc 
scalar field (indeed, it is well known that the low-energy 
limit of the bosonic string theory is an uq = — 1 Brans- 
Dicke theory [5]). Scalar-tensor cosmology, in which the 
effective gravitational coupling G g depends on time, has 
been the subject of much work [6, 7] but much less at- 
tention has been devoted to inhomogeneous solutions in 
which G e ff depends also on space. However, there is re- 
ally no support for assuming that this spatial dependence 
can be neglected [8, 9]. Spherically symmetric inhomo- 
geneous solutions of scalar-tensor gravity representing a 
central condensation embedded in a cosmological back- 
ground have been found in [9]. 

There is plenty of additional motivation for study- 
ing analytical solutions of gravitational theories repre- 
senting a central object in a cosmological space. First, 
the present acceleration of the cosmic expansion [10] re- 
quires, if one is to remain within the boundaries of gen- 
eral relativity, that approximately 73% of the energy con- 
tent of the universe is in the form of exotic (pressure 
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(m) 



) dark energy [11] (see [12] for a list of 



references and [13] for a comprehensive discussion). An 
alternative to this ad hoc explanation is that gravity de- 
viates from general relativity at large scales. Further 
motivation for alternative gravity comes from the fact 
that virtually all theories attempting to quantize grav- 



ity produce, in the low-energy limit, not general relativ- 
ity but modifications of it containing corrections such as 
non- minimally coupled dilatons and/or higher derivative 
terms. 

These ideas have led to the introduction (or better, re- 
vival) of f{R) gravity to replace Einstein theory at large 
scales [14-17] and explain the cosmic acceleration (see 
[18, 19] for reviews and [20] for shorter introductions). 
Since the f(R) theories of interest for cosmology are de- 
signed to produce a time-varying effective cosmological 
"constant" , spherically symmetric solutions representing 
black holes or central condensations in these theories 
are expected to be asymptotically Fricdmann-Lcmaitrc- 
Robertson- Walker (FLRW), not asymptotically fiat, and 
to be dynamical. Very few such solutions are known, 
among them the inhomogeneous time-dependent solution 
of Clifton in f(R) = R 1+s gravity [21, 22]. 

Second, analytical solutions describing central objects 
in a cosmological background are of interest also in gen- 
eral relativity. The first study of this kind of solution by 
McVittie [23] is related to investigations of the problem 
of whether, and to what extent, the cosmic expansion 
affects local systems (sec [24] for a recent review). In ad- 
dition to the old (and largely overlooked) McVittie solu- 
tion [23] , which is not yet completely understood [25-28] 
relatively few other solutions with similar features have 
been reported over the years [29] . 

Third, more recent interest in cosmological condensa- 
tions in the context of general relativity arises from yet 
another attempt to explain the present cosmic accelera- 
tion without dark energy and without modifying gravity. 
This is the idea that the backreaction of inhomogeneities 
on the cosmic dynamics is sufficient to produce the ob- 
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served acceleration [30]. However, the formalism imple- 
menting this idea is plagued by formal problems and it 
has not been shown to be able to explain convincingly 
the cosmic acceleration. Indeed, even the sign of the 
backreaction terms in the equation giving the averaged 
acceleration has not been shown to be the correct one 
[31-33] and recent work casts even more serious doubts 
on this proposed solution to the cosmic acceleration prob- 
lem [34]. 

Attempts to move beyond these riddles involve the 
consideration of analytical solutions of Einstein theory 
describing cosmological inhomogeneities and including 
Lemaitre-Tolman-Bondi, Swiss-cheese, and other mod- 
els [35, 36]. Moreover, the teleological nature of the 
event horizon has prompted the consideration of appar- 
ent, trapping, isolated, dynamical, and slowly evolving 
horizons ([37, 38] and references therein), a subject of 
great interest [39]. There has also been interest in dy- 
namical black hole horizons in relation to the accretion 
of dark energy [40]. Physically, black hole event hori- 
zons can only be traversed from the outside to the inside 
while, for the traditional cosmological event and particle 
horizons, signals can cross from the inside to the outside 
but not vice-versa. Event horizons are null surfaces and 
are appropriate to describe stationary situations but, as 
said, they require the knowledge of the entire spacetime 
manifold to even be defined. An apparent horizon is a 
spacelike or timelike surface defined as the closure of a 
3-surface which is foliated by marginal surfaces (those 
on which the expansion of the congruence of radial null 
geodesies vanishes) [41]. 

With all these motivations in mind, it is interesting to 
further explore analytical solutions of alternative grav- 
ity theories representing spherical objects in cosmological 
backgrounds. Here we consider the class of solutions dis- 
covered by Clifton, Mota, and Barrow [9] in Brans-Dicke 
theory, described by the action [2] 



2nL {m) 



< ^ (1) 

where n = 8itG 7 G is Newton's constant, C^ m > is the 
matter Lagrangian, and the Brans-Dicke scalar field (f> 
corresponds to the inverse of the gravitational coupling 
Geff. 1 Matter is assumed to be a perfect fluid with en- 
ergy density p^ m \ pressure p( m \ and equation of state 
p( m ) = — l) p( m ), where 7 is a constant [9]. In the fol- 
lowing sections we analyze and discuss the structure of 
the solutions of [9] , focussing on the dynamical behaviour 
of their apparent horizons, in an attempt to understand 
if the these solution harbor black holes or naked singu- 
larities. The bizarre behaviour of the apparent horizons 
we find seems to be rather typical of solutions describing 
cosmological black holes [43] in a certain region of the pa- 



rameter space, but other behaviours appear for different 
combinations of the parameters. 

Note that, even though it is standard procedure to 
rely on apparent horizons as proxies for event horizons 
to characterize black holes in theoretical and numerical 
relativity [38, 45], it is also well known that apparent 
horizons depend on the spacetime slicing adopted [44] 
(this problem is perhaps less worrisome when spherical 
symmetry is assumed) . We adopt the same practice here, 
bearing the caveat just mentioned in mind. 



II. CLIFTON-MOTA-BARROW SOLUTIONS 

We begin with the Clifton-Mota-Barrow spherically 
symmetric and time-dependent metric [9] 

ds 2 = -e^dt 2 + a 2 (t)e^ {dg 2 + g 2 dQ 2 ) , (2) 

where dVL 2 = d0 2 + sin 2 9 dip 2 denotes the line element on 
the unit 2-sphere, 
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p( m ) is the energy density of the cosmic fluid, ujq is the 
Brans-Dicke parameter, m is a mass parameter, a, </>o, ao, 
Pg 1 ™' and to are positive constants (where <f>o, /Qq™^ and 
to are not actually fully independent). Moreover, g is 
the isotropic radius related to the Schwarzschild radial 
coordinate f by 



(9) 



so that 



f = e [i + — 

2ag 



df = 1 



Aa 2 g 2 



dg . 



(10) 



1 We follow the notations and conventions of Rcf. [421 . 



The quantity a is real for uj < — 2 and for ujq > 
-3/2. For definiteness, we impose that luq > —3/2 and 
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f3 > 0. The Clifton-Mota-Barrow metric (2) is separable 
and reduces to the spatially flat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric in the limit m —> 
in which the central mass disappears. For 7^2, setting 
uj Q = (7 — 2) _1 yields j3 = and the metric becomes 
static, whereas the scalar field remains time-dependent. 
Setting 7 = 2 or 7 = 4/3 leads to (3 = 1/2 and the 
scale factor scales as y/i independent of the value of the 
Brans-Dicke coupling ujq. We will consider the physically 
interesting special cases in a separate section below. 

Our main concern here is whether the solutions in 
the Clifton-Mota-Barrow class represent black holes or 
naked singularities, embedded in a cosmological back- 
ground (by naked singularity we simply mean a timelike 
or null singularity which is not covered by an event hori- 
zon). To answer this question, we would like to determine 
the location and the nature of horizons. Since the space- 
time is dynamical, it is appropriate to consider apparent, 
instead of event, horizons and, therefore, we will locate 
the apparent horizons and study their dynamics. 



III. FINDING THE APPARENT HORIZONS 

We proceed by rewriting the metric in the more famil- 
iar form 

ds 2 = -A 2a dt 2 + a 2 {t)A^ a2 - 2 ^dr 2 + rW , (11) 
using the areal radius 

\ 2 

A ± {a -l)( a +2) 



= a(t)rAH a - 1 K a+ V . (12) 
The differential dr is related to df by 



ds 2 



A 



2a _ d 2 (t)P 2 2 (a 2 +2a _ 2) 

B(g) 2 



F 2 dt 2 
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+ ^K a _ 1)(q + 2)A i. (a - 1)(a+ 2)-2 df ) (13) 
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which means that 



df 



a{t)Ai( a ^ a + 2 ^ 2 [A 2 + ^f{a- l)(a + 2)] 



(14) 



The line element can now be written as 



ds' 



A 



2a _ d 2 (t)r 2 ,SL( a 2 + 2 a -2) 



B 2 {g) 



„ a(t)f . c 2 +3a-2 

-2-^— A — 5 drdt 



B*(q) 
B\g) 

where we have defined the positive function 



B{e) = A 2 { e ) + 



{a - l)(a + 2) m 



(15) 



(16) 



(Note that B > is a consequence of a = J > 1-) 

We now introduce a new time coordinate t which serves 
the purpose of eliminating the time-radius cross term. 
Define i such that 



dt = 



1 



F(t,r) 



[dt + i/j(t,r)dr] , 



(17) 



where ip(t, r) is a function to be fixed later and F(t, r) is 
an integrating factor which guarantees that dt is an exact 
differential and satisfying the equation 



-(- 

dr\F 



at (f 
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The line element then assumes the form 
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2a _ a 2 (t)f 2 2 (a 2 +2a _ 2) 

B(g) 2 
ipa(t)r 



„ Fd(t)f , o 2 + 3 Q -2 l 

2 - , L A — 5 \ drdt 
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The choice 



for the function ip, with 



i{t)f A- 



B 2 D{t,g) 



(20) 



D(t, g) = 1 



d 2 {t)f 2 
B 2 



(21) 



4 



turns the metric into the simple form 

/ H 2 



ds z 



-A Za DF 2 dt 2 



A 5 



is best seen using eq. (24) as at this limit A, B 1 (the 
limit is less straightforward in eq. (26) as r — > oo and 
q — > oo). Therefore, we expect the horizon at larger radii 
to be a cosmological one. 



r 2 dn 2 



(22) 



where H = a(t)/a(t) denotes the Hubble parameter of 
the background FLRW universe. We are now able to 
locate the apparent horizons (when they exist), which 
are the loci of spacetime points satisfying V c rV c r = 0, 
or g rr = [38, 46], that is 



B 4 D 



H 2 r 2 A 2 (- 2 ~ a ) +A 2 B 2 D 



0. 



(23) 



The solution of this equation reduces to the condition 
D = 0, or 



£? 2j4 2(a-i) = H-r 



2„2 



which, explicitly, reads 



A' 



A 



(a — l)(a + 2) ma(t) (a-ixa+2) 
5 A 



(24) 



= ±Hr . 



. (25) 

In an expanding universe with H > the quantity in 
square brackets is positive, hence we choose the positive 
sign. Eq. (25) can then be written as 



Hr 2 (a -!)(« + 2) 



2(a-l)(a + l) 



ma(t)A ~* ~° '"' - A Q+1 r = . 
a" 

(26) 

The Ricci scalar becomes singular as r — > for all pos- 
itive values of the mass parameter m (see Appendix A) 
and this limit denotes a central singularity. The energy 
density (8) of the cosmic fluid also diverges in this limit. 



B. The static limit 

We now consider the limit in which the metric becomes 
static, which corresponds to f3 = and yields a(t) = 
ao, see eq. (5). This value for /3 is obtained for uq = 
(7 — 2) 1 (with 7 ^ 2). This requirement implies that for 
each theory in the Brans-Dicke class, (i.e., for each value 
of luq) there is at most one solution with a static metric 
in the Clifton-Mota-Barrow family, and it corresponds 
to a specific choice of equation of state for the fluid. As 
mentioned earlier, for a to be real, one needs to have 
luq < —2 or luq > —3/2. This translates to 7 > 3/2 or 
7 < 4/3 when the (5 — condition has been imposed. 

Eqs. (6) and (8) yield 



(f>(t,r) = (f)Q 



(m) 



-2a 



(27) 



(28) 



The Brans-Dicke field <p depends on time even though 
the metric g^ v and the matter energy density p^ 1 do not. 
In fact, there is no solution in the Clifton-Mota-Barrow 
class which is genuinely static. 

In terms of the areal radius r, it is 



A 2a dt 2 



A 2 

— dr 2 
B 2 



r 2 dQ 2 



(29) 



and the apparent horizons are located by the equation 



9' 



equivalent to B = 0, or 



IV. SPECIAL CASES AND LIMITS 

Before determining the generic behaviour of apparent 
horizons, it is useful to look into some special limits, of 
cither the theory or the solutions, which will help us gain 
some intuition. 



A. The zero mass limit 

In the limit m — > in which there is no central object, 
eq. (26) reduces to Hr 2 = r, which yields r = H^ 1 , the 
Hubble horizon. 2 This value is also obtained in the limit 
of large g in which r becomes a comoving radius and the 
metric approaches the spatially flat FLRW metric. This 



-5 (a 2 -2 )q+—; =0. 



a 



Aa 2 



(30) 



The discriminant of this quadratic equation is A (a 2 ) = 

and one easily finds that A > for 

a < 1 and a > 2 (remember that a > 0, cf. eq. (7)). 
Therefore, for 1 < a < 2 there are no real roots and no 
apparent horizons. For a < 1 and for a > 2 the real 
roots, 



Q± 



m 
772 



2 ± 



(31) 



are both negative and do not correspond to apparent 
horizons. We conclude that the solution with static met- 
ric always describes a naked singularity. 



C. The limit to general relativity 



2 In a FLRW universe with curvature index k ^ the cosmological 



apparent horizon has radius (H 2 + k/a? 



We now consider the limit to general relativity ob- 
tained for ujq — > 00. When 7 7^ and 7 7^ 2, this limit 
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yields a 1, 
ds 2 = 



d 3 



00, an ' 



11 + —, 

(dp 2 + g 2 dSl 2 ) , 



2o 



(m) / *0 



(32) 

(33) 
(34) 



This metric corresponds to one of the generalized 
McVittie metrics studied in Refs. [47-49] which, in 
isotropic coordinates, assume the form 



ds 1 



1 - 



2ga(t) 



M(t) 
2ga(t) 



dt 2 



-a 2 {t) 



M(t) 
2ga(t) 



(dg 2 + g 2 dtt 2 ) , (35) 



where M(t) is an arbitrary positive regular function of 
time. 

The McVittie solution of general relativity originally 
introduced to study the effect of the cosmological expan- 
sion on local systems [23] is obtained for M(t) = const. 
This time independence of the function M(t) in this case 
follows from the McVittie condition G\ = which cor- 
responds to zero radial energy flow Tq = 0. Lifting this 
restriction and allowing for radial accretion of energy gen- 
erates the solutions (35) with general functions M{t) (see 
the discussion in Ref. [47]). It is shown in Ref. [49] that, 
in the class of generalized McVittie solutions (35), the 
solution with "comoving mass function" M(t) = Moa(t) 
(where Mo is a constant) is a late-time attractor for solu- 
tions characterized by a background universe which keeps 
expanding in the future (a — > +00). This is precisely the 
ujo — > 00 limit of the scalar-tensor solution (2)-(8), which 
also makes it clear that the Clifton-Mota-Barrow solu- 
tions are indeed accreting. Incidentally, the generalized 
McVittie solutions (35) of general relativity were derived 
two years after the discovery of the Clifton-Mota-Barrow 
solution (2)- (8) and the one with late-time attractor be- 
haviour and with M — Mq a(t) could, in principle, have 
been discovered by taking the limit to general relativity 
of this Brans-Dicke solution. 




FIG. 1: Radii of the apparent horizons in units of 
(ma,) 1 " 1- ^ as functions of time in the same units for 
ojo = —17/12. The red, dashed curve corresponds to dust 
(7 = I) and the green, solid curve corresponds to both radia- 
tion (7 = 4/3) and stiff matter (7 = 2). For dust, there is only 
one apparent horizon whose radius reaches a maximum and 
then decreases. For radiation and stiff matter, instead, there 
is a naked singularity in a universe which expands forever. 



The apparent horizons of the generalized McVittie 
metrics (35) have been discussed in [48]. For large val- 
ues of wo, the solution (2)- (8) approaches the attractor 
McVittie solution and its apparent horizons should also 
approach those of the attractor McVittie metric: jump- 
ing ahead slightly, this is indeed the case, as can be seen 
by comparing our Fig. 4 with Fig. 3 of Ref. [48] . 

The 7 = case, which corresponds to a cosmolog- 
ical constant, leads to a diverging exponent /3 for the 
scale factor a(t) when ujq — > 00. This behaviour can 
be attributed to the fact that the Clifton-Mota-Barrow 
solution assumes a power law form for the scale factor, 
whereas the general relativity limit of the solution is ac- 
tually expected to be Schwarzschild-dc Sitter spacctime. 

For 7 = 2 the ujq — > 00 limit yields a — > 1, <j) oc t~~ , 
a(t) cx \fi, p( m \t) oc t~ 3 A~ 2 , and the metric is the same 
as in eq. (32). 



GENERIC BEHAVIOUR OF APPARENT 
HORIZONS 



3 When 7 = 2 the scale factor is forced to be a(t) oc \/t, <j> oc t~ 1 
and p( m ' oc t~ :i and do not depend on the Brans-Dicke coupling 
parameter u>q. However, the limit luq — > 00 still yields a = 1 
and leads to the same functional dependence on q for the various 
quantities as the 7^2 case. The metric still belongs to the 
generalized McVittie class. 



Having discussed the special cases, we now turn our at- 
tention to the behaviour of apparent horizons in generic 
solutions of the Clifton-Mota-Barrow family. In order 
to solve eq. (26) and determine the location of these 
horizons, it is convenient to introduce the new quantity 
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FIG. 2: Same as Fig. 1 but for ojq = -1/3. The blue, dotted 
line corresponds to a cosmological constant (7 = 0). In all 
cases there is one, ever expanding horizon, and so the solution 
appears to represents a naked singularity in an expanding 



universe. 



x = in terms of which it is 

Zag 5 



A = 



1-x 
1 + x 



(36) 



while H = (j/t. One can now express parametrically the 
radius r of the apparent horizon(s) and the time coordi- 
nate as functions of the parameter x, obtaining 



r(x) 
t{x) 



fj m (1 + x) 2 ( 1 — x 
2a x \ 1 + x 



2a 



(37) 



ma*l3 (l + x)^( a+1 ) 



(1-x) 



2/a 



2x (a-l)(a + 2) {1 _ x) _ 2(a _ 1)/c 



a 



(38) 



The radii of the apparent horizons as functions of time 
are plotted in Figs. 1 to 4 for the values of the Brans- 
Dicke parameter wo = — 17/12, —1/3, 1, and 10 5 , respec- 
tively, and for various choices of the equation of state 
parameter 7. In these plots r and t are actually mea- 
sured in units of 



(ma*)' 



m \ 

a o — ) t . 

to. 



(39) 



as this convenient normalization completely absorbs the 
dependence on the parameters m, ao, to. 
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FIG. 3: Same as previously but for ljo = 1. For all three 
values of 7, at early times there is only one horizon. As the 
universe expands, the singularity gets covered by two more 
apparent horizons. Two of the horizons eventually merge and 
disappear, leaving behind only the cosmological horizon cov- 
ering a naked singularity. 



The blue, dotted curves correspond to a cosmologi- 
cal constant (7 = 0) and the red, dashed curves corre- 
spond to dust (7 = 1). The green, solid curves show 
the behaviour of the apparent horizons for both radia- 
tion (7 = 4/3) and stiff matter (7 = 2). This is because 
(3, which determines the scaling of the scale factor with 
time, is equal to 1/2 and independent of ujq for both of 
these values. For ujq = —17/12 (Fig. 1) we do not con- 
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(b) uiq = 10 , zoom-in 

FIG. 4: Same as previously but for luo = 10 5 , uiq — > oo being 
the limit to general relativity. Here for all three cases there 
are two horizons, presumably a black hole horizon and a cos- 
mological horizon. For a cosmological constant (blue, dotted 
curve) these two horizon merge and disappear. For the other 
two cases, there is initially a naked singularity which eventu- 
ally gets covered by the two horizons. 

sider the case of a cosmological constant, corresponding 
to 7 = 0, as it leads to a contracting universe. 

As can be seen in the figures (see captions for more 
details), for uj = —17/12 and wo = —1/3 there is only 
one apparent horizon for all of the values of 7 we have 
considered. In most cases, this horizon is expanding for- 
ever, so the solution is most likely to represent a naked 



singularity in an expanding universe. For ujq = —17/12 
and for dust (7 = 1), on the other hand, the apparent 
horizon exhibits a perhaps more remarkable behaviour: 
it initially expands, to reach a maximum radius and then 
contracts to reach zero radius asymptotically. 

Even more noteworthy is the behaviour of the appar- 
ent horizons when loq = 1 (Fig. 3). For dust, radiation, 
and stiff matter there is initially one expanding appar- 
ent horizon, see Fig. 3(a). Two more apparent horizons 
appear. The outer one expands, while the inner one even- 
tually merges with the initial one and they both disap- 
pear. Similar phenomenology was reported in Ref. [22] 
for Clifton's solution [21] of metric f(R) = R 1+s gravity. 4 

In fact, this puzzling behaviour was found long ago 
in the Husain-Martinez-Nuhez solution [43] describing a 
black hole embedded in a universe filled with a free mass- 
less scalar field minimally coupled to gravity and accret- 
ing onto the black hole (compare Fig. 3(a) with Fig. 1 of 
Ref. [43]). 

For ujq = 1 and 7 = 0, which corresponds to a cosmo- 
logical constant and is presented in Fig. 3(b), the situa- 
tion is similar, except for the fact that the pair of hori- 
zons actually appears inside the initial horizon. Such 
behaviour has not been reported before to the best of 
our knowledge. 

Finally, Fig. 4 corresponds to the large value of the 
Brans-Dicke parameter u>q = 10 5 . The behaviour of 
the apparent horizon dynamics is very similar to that 
present in the general relativity limit of the Clifton- Mota- 
Barrow solution obtained for luq — > 00 and discussed in 
Sec. IV C. For dust, radiation and stiff matter, the sin- 
gularity is initially naked and eventually gets covered by 
two expanding horizons, see Fig. 4(a). For a cosmolog- 
ical constant this picture is reversed: there are initially 
two nested horizons, one expanding and one contract- 
ing, which eventually merge and disappear, leaving the 
singularity naked, see Fig. 4(b). 



VI. DISCUSSION AND CONCLUSIONS 

There are relatively few solutions describing central 
matter configurations embedded in FLRW backgrounds 
in general relativity, and even fewer in alternative theo- 
ries of gravity. We have studied here the Clifton-Mota- 
Barrow class of spacetimes, which are solutions of Brans- 
Dicke theory. The latter is perhaps the minimal imple- 
mentation of a varying gravitational coupling, containing 
only a scalar extra degree of freedom. As such, it is justly 
regarded as the prototypical alternative to Einstein's the- 
ory. It is, therefore, quite interesting to assess whether 
or under which conditions can the Clifton-Mota-Barrow 



4 This fact is not surprising since metric f(R) gravity is equivalent 
to a Brans-Dicke theory with (f> = f (R), ui = and a scalar field 
potential V(<j>) [18]. 
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spacctimes describe a realistic localized matter configu- 
ration embedded in an evolving universe. 

Given that these spacctimes contain singularities, we 
have focussed our study on the behaviour of dynamical 
apparent horizons. According to the position in parame- 
ter space, we have uncovered different types of behaviour 
for these horizons. The most important result is perhaps 
that, for certain values of the parameters, the Clifton- 
Mota-Barrow spacetime appears to contain a naked sin- 
gularity (at least as far as one can tell based on the pres- 
ence/absence of apparent horizons; though unlikely, it 
is possible that the particular slicing of the spacetime 
leads to the absence of an apparent horizon even though 
the singularity is cloaked by an event horizon). In some 
cases, this singularity is present from the time of the big 
bang, thus preventing us from obtaining the metric and 
scalar field as regular developments of Cauchy data, and 
later gets covered by black hole and cosmological hori- 
zons. For other values of the parameters, pairs of black 
hole and cosmological horizons appear and bifurcate, or 
merge and disappear, a phenomenology known from a 
solution of general relativity [43] and one of f(R) grav- 
ity [21, 22]. Overall, the Clifton-Mota-Barrow class of 
solutions exhibits a great richness of behaviours of its 
apparent horizons, including the new ones reported in 
Fig. I and Fig. 3. 

The physical relevance of spacetimes harbouring naked 
singularities is, of course, questionable. However, there 
are still two scenarios in which the Clifton-Mota-Barrow 
spacctimes might still be physically relevant: (i) in the 
region of the parameter space where a black hole horizon 
eventually cloaks the singularity, it is conceivable that 
they can (approximately) describe the late time evolu- 
tion of black holes that have formed from collapse in 
FLRW spacetime (a different solution would be needed 
to describe this collapse); (ii) even in the region where 
no horizon forms, they might be able to (approximately) 
describe the exterior of a matter configuration embedded 
in an FLRW universe (a different solution will be needed 
in order to describe the interior). Whether or not any 
of these two scenarios are meaningful requires further in- 
vestigation. 

The fact that such a variety of behaviours (cos- 



mological black holes, naked singularities, appear- 
ing/bifurcating and merging/disappearing pairs of ap- 
parent horizons) is contained in the relatively simple 
Brans-Dicke theory leads us to believe that more com- 
plicated theories of gravity will exhibit an even greater 
degree of richness and complication when it comes to dy- 
namical horizons, which has not yet been explored. 

Lastly, one might be tempted to consider the ther- 
modynamics of these dynamical apparent horizons, al- 
though its physical meaning is still questioned [51]. In 
any case, it should be noted that the field equations of 
Brans-Dicke theory can be recast in the form of effec- 
tive Einstein equations G = 87r(T M „ + T^t ) in which 
the Brans-Dicke scalar field plays the role of an effec- 
tive stress-energy component T$ . The latter can easily 
violate all of the energy conditions because it contains 
terms linear in the second derivatives of 4> in addition to 
the usual terms quadratic in its first derivatives. 
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Appendix A: Ricci scalar 



The expression of the Ricci scalar is 
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Since a > 1 the Ricci scalar diverges as g — > Using 
eq. (9), it is seen that this value of the isotropic radius 
corresponds to f = 2m/a and (using cq. (12)) to the 
areal radius r = 0. Therefore, r — > denotes a central 



singularity, which is a strong one in the sense of Tipler's 
classification [50] because the area of the 2-spheres orbits 
of symmetry vanishes as r — > 0: an object falling onto 
r = will be crushed to zero volume. 
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